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^SJ ■ Pion-meson-exchange nucleon-nucleon potentials are derived for two nucleons 

. in the intermediate states. The mesons we include are (i) pseudoscalar mesons: 

ps) i TT, ry, ry'; (ii) vector mesons: p^uj^cf); (iii) scalar mesons: ao(980), e(760), /o(975); 

\ and (iv) the J = contribution from the Pomeron. Strong dynamical pair 

Q ' suppression is assumed, and at the nucleon-nucleon-meson vertices Gaussian 

■ form factors are incorporated into the relativistic two-body framework using a 

dispersion representation for the pion- and meson-exchange amplitudes. The 
Fourier transformations are performed using factorization techniques for the 
energy denominators. The potentials are first calculated in the adiabatic 
^ ' approximation to all planar and crossed three-dimensional momentum-space 

Ch ■ vr-meson diagrams. Next, we calculate the 1/M corrections. 
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I. INTRODUCTION 



This paper and its companion paper [|I| are part of our program to extend the Nijmegen 
soft-core one-boson-exchange potential [0 by including two-meson-exchange potentials. The 
material presented here is incorporated into a new extended soft-core nucleon-nucleon model, 
hereafter referred to as the ESC potential. We are still investigating ways to improve this 
model as much as possible, but the preliminary version of the ESC potential employed in the 
present paper is here used to illustrate the contributions of the various two-meson exchanges. 

In low-energy nucleon-nucleon models, exchanges with an effective mass less than 1 GeV 
seem to be the most important ones. The very heavy mesons possibly only play some role 
for S waves, but they can presumably be covered up by the lower-mass mesons. After one- 
meson exchange, the effects of two-meson exchanges at low energies seem worthwhile to be 
investigated and, as we will demonstrate with the ESC potential, they indeed provide a 
major improvement in the theoretical description of the nucleon-nucleon data. 

Recently, we have evaluated the soft two-pion-exchange potentials in coordinate space, 
using some new techniques 0. We noted that these techniques could also be applied to 
double-meson exchanges where the mesons have different masses. For the double-meson 
exchanges we can distinguish two different classes. The first class consists of the diagrams 
with two baryons (nucleons) in the intermediate states, which can be simply understood as 
the second-order contributions in a series expansion of multi-meson exchanges. In this paper 
we calculate these two-meson-exchange (TME) potentials, using Gaussian form factors, for 
the cases where one of the exchanged mesons is the pion. For exchanges with effective 
masses below 1 GeV these are the most important ones. An extension to the most general 
two- meson exchange (not necessarily involving at least one pion), can then be evaluated in 
a straight forward manner. Furthermore, we restrict ourselves to the case with two nucleons 
in the intermediate states. The potentials containing the isobar excitation effects explicitly 
have been considered elsewhere @J^. 

The second class consists of the diagrams where either one or both nucleons contains 
a pair vertex ("seagull" diagrams). The pair-meson diagrams can be viewed as the result 
of integrating out the heavy-meson and resonance degrees of freedom in the two-meson- 
exchange processes (see, e.g., Ref. p). It will appear that these "seagull" contributions can 
be interpreted successfully in a "dual" picture [0. The various possibilities for meson pairs 
coupling to the nucleon are inspired by the chiral-invariant phenomenological Lagrangians 
that have appeared in the literature. The pair-meson potentials belonging to this second 
class will be presented in the companion paper 

The mesons we include are (i) pseudoscalar mesons: vr, r],!]'; (ii) vector mesons: p,uj,(l); 
(ii) scalar mesons [^: ao(980), £:(760), /o(975); and (iv) the J = contribution from the 
Pomeron. The 27r-exchange potential has already been derived in Ref. 0, and we here 
only include those results for reasons of completeness. In the literature, the treatment of 
two-meson exchanges is incomplete. For example, the iruj potential has been derived in 
the dispersion-relation approach [Q, which does not allow for the inclusion of Gaussian 
form factors as applied in the present study. For the ttuj, vrp, and vrcr exchanges, only 
the non-iterative potentials from the graphs with at least one meson in the intermediate 
state (see below) have been evaluated in momentum space in the framework of noncovariant 
perturbation theory 0]. In this paper we give, apart from the 27r-exchange, the potentials 
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for all pion-meson exchanges with a total mass less than ~ 1 GeV. 

Our general approach to two-meson-exchange potentials is given, in principle, in Ref. [Q. 
Starting from the relativistic two-body equations, the two-meson-exchange potentials for 
the relativistic three-dimensional integral equation and for the Lippmann-Schwinger equa- 
tion are derived from the second-order Feynman diagrams by integration over the relative 



energies |]rT|. This way, the "old-fashioned" perturbation graphs for the TME potentials are 
obtained from the Feynman diagrams. 

The proper way to incorporate Gaussian form factors is explained in and, as al- 
ready mentioned above, in this paper the important "factorization" technique involving two 
mesons with different masses will be exploited. It is in particular this factorization technique 
which enables us to express the two-meson exchanges in terms of the basic coordinate-space 
function 0c'(^) Nijmegen one-boson-exchange (OBE) model P JT^JT3[] . Technically, 

there is no compelling reason to work in coordinate space, rather than in momentum space. 
Electromagnetic effects can be included more easily in coordinate space, whereas relativis- 
tic effects can be studied more effectively in momentum space. However, we think that 
at low energies relativistic effects can be included sufficiently accurate in coordinate space. 
Coordinate-space interactions are useful, for example, in Jastrow variational calculations 
(see, e.g., Refs. |T^JT^) or in applications where one wishes to distinguish a (long-range) 



meson-exchange region and a (short-range) quark-exchange region (see, e.g., Refs. [16,1] 
Furthermore, the effects of the particular shape of the central, spin-spin, tensor, and spin- 
orbit components on the phase shifts is much more transparent in a coordinate-space picture. 
Ultimately, it is our goal to produce an NN potential which is exactly equivalent in both 
coordinate space and momentum space. 

The graphs which we calculate are: (i) the planar and crossed vr-meson-exchange graphs 
of the similar type as those calculated by Brueckner and Watson for two-pion exchange, 
and (ii) the iterated vr-meson-exchange graphs of the type of graphs calculated by Taketani, 
Machida, and Ohnuma for two-pion exchange. As the distinction between these two 
classes of graphs is convenient as a means to denote the different contributions, we will also 
adopt in this paper the following nomenclature: the non-iterative potential graphs with at 
least one meson in the intermediate state are of class (i) and will be referred to as planar 
(/ /) and crossed (X) BW graphs, while the iterative graphs with only nucleons in the 
intermediate state are of class (ii) and will be referred to as TMO graphs. 

The vr-meson-exchange potentials are calculated first using the adiabatic approximation 
for the nucleons in the intermediate states. However, in the meson exchanges we include the 
1/M^ terms, in analogy with the OBE model p|,|12|,|13[. We find that up to this order the vre 



and vr/o potentials are identically zero. Also, in the vr-vector and vr-pseudoscalar potentials 
there are large cancelations between the contributions from the planar and crossed-box 
diagrams when the meson is an isoscalar (J = 0) meson. Therefore, the next to leading 
order is expected to be important, and so we here also calculate the nonadiabatic 1/M 
corrections and the 1/M corrections from the pseudovector coupling of the pion and the 
other pseudoscalar mesons. 

The paper is organized as follows. In Sec. || we define the two-meson-exchange kernels. 
We give the interaction Hamiltonians and briefly indicate how the Gaussian form factors can 
be implemented. In Sec. |II|the definition of the nucleon-nucleon two- meson potential for the 
Lippmann-Schwinger equation is briefly repeated from |3] and the vertices in Pauli-spinor 
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space for the different couplings are given. Here we also mention the approximations made 
in these vertices. In the next section, the vr-pseudoscalar, the vr-scalar, the 7r-Pomeron, and 
the TT- vector potentials are derived. In Sec. [V| we derive the 1/M corrections, due to (i) the 
nonadiabatic corrections, (ii) the 1/M term in the pseudovector coupling of the pion and 
the other pseudoscalar mesons, and (iii) a cancelation between the off-shell term in 

the pseudovector vertices and the 1/M term of the two-nucleon intermediate state in the 
TMO diagrams. In Sec. |V| the results are shown and discussed. 

Appendix ^ contains a dictionary of differentiation formulas, adequate for deriving the 
final form of the potentials of this paper in coordinate space. Explicit expressions for the 
leading-order contributions to the coordinate-space potentials are given in Appendix O. 



II. THE TWO-MESON-EXCHANGE KERNEL 

The fourth-order two-meson-exchange kernel is derived following the procedure as dis- 
cussed in , to which we refer for details and definitions. The only difference is that here we 
choose to shift the M/ E factors into the Green's function. This means that the Thompson 
equation reads 



(p') = 0?1(p') + 4^Hp';W^) 

X |f/3pir-(p',p|W^)0++(p) 



where now the Green's function is 



(2vr) 



MaMk 



E,{^')E,{^') 



a^(pOa^(-pO 



(2.1) 



(2.2) 



and so the irreducible kernel is then given by 

ir-'^(p', v\W) = -{27ir^[W - W(p')] [W - W(p)] / dp', / dpo 



X 



X 



-1 



-1 



(2.3) 



The transition from the Thompson equation ( p.l| ) to the Lippmann-Schwinger equation, 

0(p') =0(°Hp')+^(p';W^) ld'pV{p',p\W)<Pip), (2.4) 
is effectuated by the transformations 

0++(p) = A^(p;iy)0(p), 

X'"(p', p\W) = N-\p'; W)V{p', p|l^)Ar-^(p; W), 

Et\p; W) = N^{p; W)gip; W), (2.5) 



with the nonrelativistic Green's function 
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where p, is the on-shell momentum. For AW scattering the equal-mass approximation is 
very nearly valid. In that case we have, taking the nucleons on the energy shell, 



N{^-W) = ^M/E{^). (2.7) 

Taking further into account the normalization factors of the nucleon Dirac spinors, it is then 
natural to write 



np',p) 





(E{^') + M\ 




\ 2M ) 



^ IE{p) + M\ I M 



2M ; V ^(p) 

^U{p',p)[l + 0{1/M% (2.8) 

where the potential function f/(p',p) has a reduced energy dependence. 

The two-meson-exchange kernel is written as a power series in A, which denotes the 
number of nucleon-nucleon-meson {NNm) vertices. The various two-meson contributions 
to the fourth-order A terms of the two-meson-exchange (TME) kernel define the various 
two-meson-exchange potentials. The corresponding fourth-order elastic AW matrix element 
of the kernel now reads 

K('\p', p\WU';ab = -{2n)-'[W - W(p')] [W - W(p)] 

^ X! J^P'o J^Po Jdkio Jdk2o j dki J dk2 i{2Ti)~^5^{p' - p - ki - 

a" ,b" 

[[T,F^\p + ki,po + k,^)T,]^^"\T,F^\-p - kl, -po - k,^)T,f'^ 

+ [T,F^\p -f ki,po + fcio)r.]('^"Hr, F^\-p' + kl, -p'^ + A:io) r,](^")} 



X 



X 



fI^\p,Po)F^\-P, -Po)V [kl -ml + i5r\ (2.9) 



Here mi and m2 denote the two meson masses, and the Fj and Tj denote the nucleon- 
nucleon-meson vertices, which follow from the interaction Hamiltonians (see below). Because 
here we only consider nucleons in the intermediate state, we have a = a' = a" = N and 
h = h' = h" = N . Note that the first term between the curly brackets corresponds to the 
planar-box two-meson-exchange diagram and the second term to the crossed-box two-meson- 
exchange diagram. In these diagrams only the contribution of the positive-energy nucleon 
states are included, in accordance with the pair-suppression hypothesis that we used in our 
earlier work on two- meson exchange 

The fourth-order potential V^'^'> consists of two parts. The first part is given by the fourth- 
order planar and crossed BW diagrams, shown in Figs. [I|(a) and [I|(b)-(d), respectively, and 
their "mirror" counterparts. The second part is due to the fact that we do not distinguish the 
two different time orderings (i.e., the explicit energy dependence) in the one-boson-exchange 
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diagrams. Hence, solving the scattering equation does not generate the time-ordered TMO 
diagrams of Figs. |^(a) and @(b), nor their "mirror" counterparts. We therefore include these 
TMO diagrams explicitly |20|, subtracting the once- iterated one-meson contribution 



K^L = K^'^ = V^'^g V^'\ (2.10) 

and so the fourth-order potential reads 

V{p', p) = K^^^ - K^^^g K^^\ (2.11) 

The procedure to derive the kernels for the planar and crossed graphs corresponding 
to Figs. ^ and ^ is amply described in Refs. [^|[ and will not be repeated here. From 
the worked-out vertices (see below) and closely following the procedure as discussed in the 
referred papers, the contribution of the different graphs can be written down immediately. 
Finally, when the mesons are distinguishable, we also have to include the contributions from 
the diagrams where we interchange the two meson lines (dashed and dotted lines in Figs. |l| 
and They merely give rise to a factor of two when we evaluate the potentials. 

For point couplings the nucleon-nucleon-meson Hamiltonians are |21[ 



npv = ^^iJ-i^-i^T^-d^'ciypv, (2.12a) 
ns = gs^rij-<t)s, (2.12b) 
T<v = gv^li,'r^-4>v - ^^^t^y'^^-^'^^v^ (2.12c) 

where (f) denotes the pseudoscalar-, scalar-, and vector-meson field, respectively. For the 
isoscalar (/ = 0) mesons, the isospin Pauli matrices, r, are absent. The nucleon-nucleon- 
Pomeron vertex has the same Lorentz structure as the isoscalar scalar vertex. 

As is well known, the pseudoscalar one-meson exchange can also be represented by the 
pseudoscalar interaction Hamiltonian 

Hps = gpsfpijbri) ■ (fips- (2.13) 

Both Tipv and "Hps give rise to the same on-energy- shell one-meson-exchange potential 
provided that the coupling constants satisfy the relation /py/m^ = gpsjlM. From this we 
tacitly assume that the pseudovector coupling constant jpy is C(l/M). This simphfies the 
expressions for the pion-pseudoscalar potentials considerably. 

The generalization of the interaction kernels to the case with a Gaussian (or any other) 
form factor has been treated and explained in 0. The same procedure can be applied to 
the various meson exchanges to be considered here. The form factors FJ^^ and F^(k^), 
which describe the meson-exchange and vr-exchange amplitudes, respectively, are simply a 
product of the Gaussian vertex form factors F^^a^)^^ = exp(— k^/2A^^^) and F^Nni}^'^) = 
exp(-kV2A^^J. 



III. DEFINITION OF NUCLEON-NUCLEON 
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TME POTENTIAL 



The transition from Dirac spinors to Pauli spinors is reviewed in Appendix C of Ref. 
There we derived the Lippmann-Schwinger equation, 



Xip) = X^'\p) + ~9{P ) Jd'p V(p',p) xip) 



(3.1) 



for the Pauh-spinor wave functions x(p)- The wave function x(p) cind the potential V(p', p) 
in the Pauh-spinor space are defined by 



(P) = XaaaSP) Ua{p,(Xa)Ub{-p,ab) 



(3.2) 



and 



x^'xi:C V xS^xS = ^^(p'' OM-p', ^i)vip', p) 

X Ua{p,aa)ub{-p,ab). (3.3) 
Like in the derivation of the OBE potentials P,p!2|, we make the approximation 

E{p) = (p2 + m2)1/2 ~ m + pV2M 

everywhere in the interaction kernels of Sec. |V], which of course is fully justified for low 
energies only. We have a similar expansion for the on-shell energy, 

W = 2(p2 + m2)1/2 ~ 2M + p2/M. 

In contrast to this kind of approximation, the full dependence of the form factors is kept 
throughout the derivation of the two-meson-exchange potential. Note that the Gaussian form 
factors suppress the high-momentum transfers strongly. This means that the contribution to 
the potentials from intermediate states which are far off-energy-shell cannot be very large. 
It also means that the 1/M expansions can be expected to be valid up to larger energies as 
well. 

The reduction of the TME potential from Dirac-spinor space to Pauli-spinor space is 
completely similar to the procedures discussed in Refs. IQ-^ . The vertex operators in Pauli- 
spinor space up to order are given by 



u{p')Tfsu{p) = 
u{v')Tfu{p) = gs 
uip')T^^Kip) = gv 

gv 



'I 



PV 



o-i-k 1- 



p'^ + p^^ 

8M2 



±^-r(p+p) + ^^-i-(p+p) 



gps 

2M 



J2 



O-i-k 1 



P^ + P 

8M2 



P 



12 



P^ 



(Ti- 



[P' + P) 



1 + 



p' p + icTi p' X p 
p' p + icTi p' X p 



(3.4) 
(3.5) 

(3.6) 



Ky- 



(p' - P)^ - 2io-i-p' X p 



2M 



4M2 

p' + p) + i{l + Kv)(Ti X k 



4M2 



V-, 



(3.7) 
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where k = p' — p and = fv/dv- The expressions for u{—p')V^'^^u{—^) are trivially 
obtained by substituting (p', p, k, cti) (— p', — p, — k, cr2). For the isovector (/ = 1) 
mesons, the appropriate isospin matrices have to be attached. In the F-matrix element of 
Eq. (|3.4|), the upper sign applies to the creation and the lower sign to the absorption of the 
pseudoscalar meson at the vertex. 

Although we will assume the pseudovector coupling for the pion, we here have also listed 
the vertex operator in case of pseudoscalar coupling for reasons of completeness, which might 
be helpful for later work by other groups. The changes can be easily accounted for, as we 
will indicate in the following sections. 

Useful for the evaluation of the second-order diagrams are the relations 

TjTi = 6ij + iejikTk, aja, = 6ij + iejikCTk. (3.8) 

Products of this type will occur for each nucleon line. The isospin factors for the planar (/ /) 
and crossed (X) diagrams can readily be evaluated. In the case of two isospin-1 mesons, we 
have 

d£\l)=3-2TvT2, 

= 3 + 2ri-T2, (3.9) 

where / denotes the total isospin of the NN state. For the exchange of an isospin-1 with an 
isospin-0 exchange, one has 

C^^ffil) = C^NJil) = C^Nk'\l) = ri-T2. (3.10) 

In order to obtain the contributions to the potentials in the adiabatic approximation, 
we expand the energy denominators in the expressions for the planar and the crossed-box 
diagrams and keep only the leading term. So, for example. 



1 1 



E{p) + E{p") -W+LU LU 



kik2-q(ki-k2 

2Mu 



(3.11) 



where we have made the on-energy- shell approximation p — p^ = and p — p^ = 0, with 
Pi the initial-state momentum, and where we introduced q = |(p' + p)- The second term 
between square brackets gives the nonadiabatic contribution to the potentials, presented in 
Sec. |V A| . The term proportional to q-(ki — k2) usually drops out for symmetry reasons. 



The only exceptions occur for the planar-box nonadiabatic pion-pseudoscalar potentials (see 

Sec. [yil) . 

In the following, we will concentrate on the derivation of the pion-meson potentials, 
which are expected to be the most important ones. The derivation of the more general case 
of two arbitrary mesons (not necessarily involving at least one pion) will then be straight 
forward. 



IV. PION-MESON-EXCHANGE POTENTIALS 
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A. Pion-pseudoscalar exchange 



Using the assumption that the pseudovector couphng constants are 0{1/M) [see the dis- 
cussion below Eq. ( p.l2a| )], we can neglect all contributions and the pion-pseudoscalar- 
exchange potential is the simplest potential we will encounter. As a representative of the 
pseudoscalar mesons we select the rj and derive the vrr^-exchange potential. The vrvr-exchange 
potential has already been discussed in Ref. ||^. The purely off-shell term proportional to 
(p'^ — p^) in the pseudovector vertex ( |3.4D gives rise to a 1/M term in the TMO diagrams, 
rather than a term. We defer this contribution to the potential to Sec. [VC| . Similarly, 
the contribution due to the oj /M term will be discussed in Sec. [VB| . 

We assign (kijCJi) to the vr meson and (k2,ci;2) to the rj meson. The planar BW graph 
of Fig. |l|(a), the TMO graphs of Fig. ^, and their "mirror" counterparts give 



fNNr)\ ( InNtt^ 



X 



ki-k2 - ^^^(ki X k2 



7) //^|^^''^"'^'^^'^(k?)^.(k: 



ki-k2 - icr2-(ki x k2) 



D//{uJi,UJ2), 



(4.1) 



and the crossed BW graphs of Figs. |l|(b)-(d) and their "mirror" graphs give 

2 / r \ 2 



X 



ki-k2 - 'icri-(ki x k2) 



7r)6 

ki-k2 + icr2-(ki X k2 



(4.2) 



FLr^ikj) and F,(k2) 



^^^^(kl), and the energy denominators Di{uJi,uj2) 



Here F^{k^) - ^-^^^y 

are given in Table |. Choosing the pseudoscalar coupling (|3.5|) instead of the pseudovector 
coupling ( |3.4| ) will result in the same leading-order potential. 

The full separation of the ki and k2 dependence of the Fourier integrals can be achieved 
as shown in Ref. 0. Using the A-integral representation, Eq. (B.4) of 0, one can derive 
that 



and 



1 2 




■ 1 




1 ■ 






Jo A2 






+ A2J 




1 1 


2 r^dX 


" 1 


1 


Ujfu!2 UJi + ijJi 


TT Jo 


X 


' 1 


uj! + X\ 

1 








ujI + X\ 



(4.3) 



(4.4) 



In performing the Fourier transform, the term gives rise to the basic function l2{m, r) 
{m/4:7i)(j)1.{m, r) , where 



m, r 



-"^^erfc - ^) 



-e"^"erfc 



A 

m 

A 



+ 



2 / 

Ar 



2mr 



Similarly, the 1/(^7^ + A^) term gives rise to the function 



F„(A,r) = e-" h{J ml + \\ r) . 



(4.5) 



Therefore, the typical coordinate-space function for the planar and crossed diagrams is of 
the form 



2 f^dX 



Ba/3{ri,r2) = - I ^ h{ma,ri)l2{mp,r2) 
71 Jo I 

-F^{X,r,)Fp{X,r2 
Accordingly, we derive from Eqs. ( f4.1| ) and ( [4. 21 ) the form 



TTTI 



i=//,x \ 



fNNv\ ( f. 



NNtt 



lim 

ri,r2-^r 



xO«(-^Vl,-^V2) i?^L.(n,r2), 



where 



B 



(//) 

NNnriy 



(4.6) 



(4.7) 



(4.8) 



and the operators Op^. (ki, k2) are given in Table ||. We should point out that these operators 
contain the contributions from all graphs, including the contributions where we interchanged 
the TT and 1] lines. This gives an overcounting when the particles are identical, and so in case 
of TTTT exchange these operators should be divided by 2. 



B. Pion-scalar exchange 



To be definite we select from the scalar mesons the ao(980) and derive the vrao-exchange 
potential. The 7r£(760)- and the 7r/o(975)-exchange potentials can simply be obtained by an 
appropriate change of the isospin structure. In the following, (ki,ui) refers to the tt meson 
and (k2,ci;2) refers to the Qq meson. 

We will evaluate the potentials up to order Using the vertices of Eqs. ( p. 4]) and 



(3.6), and including the normalization factors of the nucleon Dirac spinors in the intermediate 
two-nucleon state, which also contribute a factor of order we readily obtain the 

contributions to the NN potentials. After some rearrangement, the planar graphs give 



V!^l{//) = Ciiff{I)g 



2 ( InNtt 



X 



p' p" + icr i p' X p" 



4M2 



o-rki)(tT2-ki)iV//(p") 



1 - 



p'-p" + icTa-p' X p" 



4M2 



8M2 



<7i-kicr2-(p" + p) + cri-(p" + p)cr2-ki 



D//{ui,U2), 



(4.9) 



and the crossed graphs give 
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1 



2 

AWao 



7. 



X 



p'p" 



AWtt 



icTi p' X p" 



4M2 



(27r)6 

>rki)(cr2-ki 



1 - 



p" -p + icTa p'" X p 

4M2 



iVx(p",p" 



+ 



_,///2 „//2 
P 1 / ' I w\ I P 



P^ 



cri-(p" + p)cr2-ki LDx(a^i,o^2). (4.10) 



8M2 ^ ^ " ' ^ ^ ' 8M2 
Here we introduced momentum vectors for the intermediate nucleon lines, given by 



p" = p + ki = p' - ks, 



P 



P 



p' - ki, 



(4.11) 



and the functions 



1 + 



P 



112 



P^ 



4M2 



iVx(p",p'") = l + 



P 



//2 



8M2 



2 ///2 

+ 



P 



12 



M2 



(4.12) 



Again, the energy denominators -Dj(co'i, ^72) are given in Table |, and for all graphs ^^(k^) 
F2^,(k?) and F,„(k2) - n.2 
coordinate-space potential can be written in the form 



-^AWao(k2)- Including the "time- reversed" diagrams, the resulting 



y(o) 



Trao 



2 

AWao 



7. 



lim 



i=//,X 

xO«(-zVi,-^V2) i?SLao(n,r2), 



AWtt 



(4.13) 



where i?^]y^„y(ri, r2) is defined by Eq. ( |4.^ ). The operators 05''(ki,k2) can be found in 
Table |I|. They contain the contributions from all graphs. The choice of pseudoscalar 
coupling ( p.l3| ) for the pion instead of pseudovector coupling ( |2.12a| ) results in a change of 
sign in one of the terms, as indicated in Table ||. 

The purely off-shell contribution proportional to (p"^ — p^) gives rise to a 1/M term in 
the TMO diagrams, which will be discussed in Sec. |V C| . Similarly, the contribution due to 
the uj/M term in the pion vertex (|3.4|) will be discussed separately in Sec. |V B. 



Note that, up to this order, both planar and crossed diagrams give rise to the same 
momentum operators Os'(ki,k2). As a consequence, the planar and crossed contributions 
for the Tie and tt/o potentials exactly cancel, and hence 



y w = V 



(0) 



0. 



(4.14) 



C. Pion-Pomeron exchange 

The pion-Pomeron potential has the same structure as the pion-scalar potential, except 
for an over-all minus sign. Hence, 



11 



Vi? = 0. (4.15) 

Although the leading-order part of the potential vanishes, the 1/M corrections to the 
potential do not, and so in order to be able to evaluate the pion-Pomeron potentials in coor- 
dinate space, it is appropriate here first to point out the peculiarities of Pomeron exchange 
with respect to meson exchange. We use as a working hypothesis that there is a contribution 
to the potential which corresponds to the Pomeron phenomenon at high energies. In QCD, 
the physical nature of Pomeron exchange is understood as color-singlet two- or multigluon 
exchange [p^ -p^. The assumption that the Pomeron for our purposes can be described 
effectively by a simple Regge pole is part of our phenomenological and practical approach. 
In the Nijmegen OBE models Pomeron exchange provides a significant contribution to the 
short-range repulsion. We refer the interested reader to Refs. P, p!2|j25[] , and references cited 



therein. It remains to be seen whether Pomeron exchange, using some plausible form for it, 
is still a useful concept when we include many more contributions (i.e., two-meson contribu- 
tions) than were taken into account in the Nijmegen OBE models. But the inclusion of the 
potentials involving the Pomeron is important in order to have the full consequences of the 
Pomeron worked out. Hence, at this stage these potentials are important for an extension 
of the Nijmegen OBE model. 



We treat the Pomeron as a Regge pole, using the Durand-van Hove model Accord- 
ingly, in dealing with Pomeron exchange, we replace the OBE factor for meson exchange by 
an infinite sum over all even J, i.e. 

E(2^+l)g|. (4.16) 

J even J \ / 

where 

ujj{t) = m'^j-t^k^ + 712]. (4.17) 



Following the same procedure as in the OBE model ||2^ (i.e., the Khuri- Jones procedure for 
a Regge pole), and using the Sommerfeld- Watson transform, one gets for the Durand-van 
Hove sum in Eq. ( [4.16| ) the expression 

Ap{t)= Y: {2J+l)Aj{t), 

J even 

^--W = -''W^fipT- (4-18) 

Here a{t) is the solution of the equation rrij — t = 0, Pit) is the residue at the pole, and ^ 
is a parameter in the Khuri- Jones representation (see [^). For low energies, our region of 
interest here, ^ is such that one can restrict oneself to the lowest J value in Eq. ( [4.18 ), which 



is J = 0. However, we should point out that in the case of the Pomeron the J = contribu- 
tion from the Khuri- Jones procedure has nothing to do with a physical particle. With the 
Chew ghost-killing factor, f]{t) = —a{t) exp(apt), one obtains for Pomeron exchange at low 
energies from Eq. ( [4.18| ) the representation 
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where the Pomeron-nucleon-nucleon couphng and the effective Pomeron mass can be ex- 
pressed as 



(0) {srIM' 



2 /^x /- ^9\"-P(0) 

9p = 1p 

mp = - ap + ttpln— — 



(4.20) 



With this preparation, we can now explain the treatment of the pion-Pomeron contribu- 
tion to the potentials. As we will see in the next section, for the Pomeron contributions we 
encounter in the energy denominators the factors l/uj'j and l/cuj. The handling of the first 
term is discussed above and that of the second is easily inferred by observing that 

^.w= = (4.21) 

J even J \ / 



From Eq. ( [4. 19|) we therefore find that 



1 



Mt) = iTTYMt)- (4.22) 



Anip 



The result of this analysis is that the basic Pomeron functions are 

/2(m,r) lGimp,r), 
1 



him, r) ^ ^/^(mp, r), (4.23) 



Amp 



where 



With the substitutions ( |4.23|) , the treatment of the Pomeron is now completely analogous 
to that of the mesons. 



D. Pion- vector exchange 



Due to the complexity of the vector-meson vertices (p.7|) , it is not very illuminating to 
give the various intermediate steps in the evaluation of the pion-vector-exchange potential. 
It suffices to note that it is convenient to separate the potential in four parts proportional to 
1, (1 + Ky), (1 + 2kv), and (1 + Ky)^, respectively. We refer to these as the electric-electric 
(e, e), the two electric- magnetic (e, m) and (e, 2m), and the magnetic- magnetic (m, m) terms. 
Assigning again (ki,u;i) to the vr meson and (k2,u;2) to the p meson, we can then write the 
resulting potential as 



i=ll,x 



NNtt 

m^ 
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oil + o(:L + o«^}(ki,k2) 



1 



2M2 

+0«(ki,k2 



(4.25) 



where the operators 0*^*''(ki, k2) can be found in Table 0. They contain the contributions 
from all graphs. Again, the choice of pseudoscalar couphng for the pion instead of pseu- 
dovector couphng involves a simple change of sign in one of the terms in O^*^. 
The coordinate-space potential reads 



NNp 



i=//,X 



NNtt 



lim 

ri,r2— >r 



X 



2M2 

+0«(-^Vi,-zV2) 



BNNnpiri,r2) 



(4.26) 



The and 7r0 potentials are again obtained by an appropriate change of the isospin struc- 
ture. Note that in leading order (no l/M^ contributions), the planar and crossed-box con- 
tributions exactly cancel for both ttuj and 7r0 potentials; see Eq. 



V. 1/M CORRECTIONS 
A. Nonadiabatic corrections 

The nonadiabatic correction from the 1/M expansion of the energy denominators is 
explained in Ref. |^. The expansion of the energy denominators involves a momentum 
dependence which for all cases can be rewritten in the form [ki-k2 ± q-(ki — k2)]/2M. 
Separating-off the momentum dependence^ the energy denominators of Table | give rise to 
the nonadiabatic energy denominators as listed in Table |T1[ Note that = ~'^^'^fj i 
and so now the contributions from the planar and crossed diagrams for 7r£:(760), 7r/o(975), 
and ttP do not cancel. We should point out that this expansion is only approximately valid, 
and in principle breaks down at the pion-production threshold (Tiab ~ 280 MeV). 

In the crossed-box diagrams the q dependence always cancels. Similarly, for pion-vector 
and pion-scalar (pion-Pomeron) exchange, the interchange of the two meson lines also re- 
moves the q dependence. Hence, for pion-scalar exchange we get 

enO^C.-(/).^..(^)\l) 

X (krk2)(o-i-ki)((T2-ki)A-'H^i,^2), (5.1) 
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where we took the vre potential as a particular example. Obviously, for pion-Pomeron 
exchange we get the same expression as for pion-scalar exchange, except for an over-all 
minus sign. For pion-vector exchange we get 



'I 



NNtt 



1 

M 



X (ki ■ ks) (o-i ■ ki) (0-2 ■ ki)Df ^ (^1, UJ2), 



(5.2) 



where we took the vrp potential as a particular example, and only kept the leading term in 

Og(ki,k2). 

The situation for pion-pseudoscalar exchange is more subtle. As mentioned above, in the 
crossed-box diagrams the q dependence cancels, resulting in 



fNNr)\ ( fNNTr\ f 1 



d^kid^k 



(27r) 



X (ki-ka 



^ki-ka)^ + cri-(ki x k2)cr2-(ki x kg) 



(5.3) 



where we took the tit] potential as a particular example. In the planar-box diagram, on the 
other hand, the q dependence survives and, after rearranging, we are left with 



^ c-s'ii) (^) (^) (^) //^e--' FM)F.(^i) 



X (ki-ka) 



'ki-ka)^ - cri-(ki x k2)cr2-(ki x ka 



+i{(Ti + (T2) ■ (ki X k2)q- (ki - k2 



(5.4) 



All integrals in Eqs. (|5rT|) - (|5l^) can be readily evaluated using the results of Appendix 0. 
Inspection of the energy denominators dI:^\uji,uj2) reveals that we need derivatives of the 
function l4^{m,r), which is defined by 



(5.5) 



while Ii{mp,r) for the Pomeron is given by the second expression in Eq. ( |4.23| ). 



B. Pseudovector-vertex corrections 

The pseudovector vertex gives rise to 1/M terms as shown in Eq. ( p.4[ ). In the planar-box 
diagrams, the 1/M contribution from each diagram is cancelled by the 1/M contribution 
from its "mirror" diagram, and therefore we only have to consider the crossed-box BW 
diagrams. The pseudovector-vertex corrections are very simple to derive for the pion-scalar 
and pion-vector potentials. Below, we again take a representative meson for each type. 
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BccaTisc of their simplicity, we also give the explicit coordinate-space expressions for these 
potentials. For pion-scalar exchange we get 

/ f \ 2 

T/(16) _ ^a:0)/.x 2 I JNNn 



X 



X 



(27r)6 



1 

M 

F^(k2)F,(k: 



((Ti-ki)(<r2-k2) + (<Ti-k2)(<r2-ki) 



_r^(1.0)/7-N„2 



'I 



NNw 



3M 



xl2{mn,r)l2{me,r) [(cti-cts) +5*12] 
For pion- vector exchange we get 



(5.6) 



1 

M 



X 



X 



(2vr) 



(cri-ki)((T2-k2) + (cri-k2)(cr2-ki) 



2a;ia;2 



3M 



xl2(m„,r)l2(mp,r) [(<ri-o-2) +5'i2]. 



(5.7) 



The pion-pseudoscalar case is more complicated than the pion-scalar and pion-vector 
cases, because now we also have to include the 1/M corrections from the pseudovector 
vertex of the 77 or rj' meson. Again, the contributions from the planar-box diagrams are 
cancelled by their "mirror" diagrams. Including the "time-reversed" diagrams, and defining 
the function 



/o(A,r) = 
we find for the crossed diagrams 



1 



1 

M 



47r2V7r 

(Pki(Pk2 i(ki+k2)-: 

(27r)6 



(5.8) 



F.(k?)F,(k^) 



x<^ (krk2)(k2 + k^) + i(o-i + 0-2) • 



qx (ki + k2)(ki-k2) + (ki x k2) q-(ki - k2) 



fNNri\ ( fNNTr\ 



1 

M 



X 



+ 



(ml + mj) l2{m„, r)/2(m^, r) - /o(A7r, r)l2{m^, r) - I'2{m^, r)/o(A^, r) 
I2{m^,r)l2{rnrj,r) + I2{rn^,r)l2{rnr,,r)^ + ^I2{m^,r)l2{mr„r) L-s|. 



(5.9) 
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C. Off-shell corrections in TMO diagrams 



Next to the type of nonadiabatic (off-sheU) corrections discussed in Sec. [V A| , there are 
1/M corrections in the TMO diagrams, due to off-shell factors like (p"^ — p^) from the 
pseudovector vertices [see Eq. (p.4|)j. In the on-shell approximation, p'^ — Pi = P^ — Pi = 0, 
this type of momentum dependence in the numerator cancels the p"^ — pf momentum 
dependence of the two-nucleon intermediate state in the energy denominator of the TMO 
diagrams, leaving a 1/M correction to the potential. 

For the pion-scalar and pion-vector exchanges it can be easily shown that these (what 
we call) off-shell TMO contributions involve the momentum operator 



Ooff(ki,k2) 



M 



CTi-kiCTa-ki =F ^(cri ■kicr2-k2 + cri-k2cr2-ki 



(5.10) 



where the =F sign refers to pseudovector or pseudoscalar coupling for the pion, respectively. 
Hence 



'{1)9 



NNe 



NNtt 



( 1 



xF^(k2)F,(k2)Ooff(ki,k 



V2M 

. 1 



' 2 2' 



(5.11) 



T/{ic) _ r<(//)|'/-^„2 



NNw 



f 1 



\2M 



xF^(k2)F,(k^)Ooff(ki,k2) 



(5.12) 



Because the off-shell factors in the pseudovector vertices occur with the same (cti, (T2) spin 
dependence as the uj/M terms, it should be obvious that the resulting potentials are very 
similar to the pseudovector-vertex corrections derived in the previous section (except for a 
relative minus sign in the spin-orbit pion-pseudoscalar part). We find 



I / fNN-K\ f 1 



+ 



\2M 

LS 



T" T' 



T' T" 



-T' T' 



VI. RESULTS AND DISCUSSION 



(5.13) 



The complete pion-meson-exchange potential can be written as a sum of all pion-meson 
exchanges ^(vra), where 



n=0, la, 16,1c 



(6.1) 
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Here, the second meson is denoted by a. Each potential consists of central, spin-spin, tensor, 
and spin-orbit parts. The leading-order contributions to the coordinate-space pion-meson 
potentials are given explicitly in Appendix 0. The nonadiabatic 1/M corrections can be 
easily derived from Sec. |V A| and Appendix 0, the pseudovector 1/M corrections are given 
in Sec. |V B| , and the extra 1/M corrections in the TMO diagrams are given in Sec. |V C] . 

In the following we show the various contributions to the potentials, using the meson- 
nucleon coupling constants and cutoff masses as given in Table They correspond to the 
values of a preliminary version of the Nijmegen extended soft-core (ESC) potential which, 
next to the standard one-boson exchanges and the two-meson exchanges discussed in this 
paper, also contains the pair-meson exchanges to be discussed in the following paper [|I|. 
In this model, the pion-nucleon coupling constant at the pion pole was fixed at /^^^/47r = 
0.0745, in accordance with the recommended value |2^. Keeping in line with the philosophy 
of the Nijmegen group p|,|T2|,|l3| , we use SU(3) relations for the coupling constants. However, 
since the existence of a scalar nonet (and hence its quark content) is still controversial, the 
three nucleon-nucleon-scalar coupling constants are fitted independently. To further limit 
the number of free parameters, the value for {f / g)p was fixed at the vector- meson dominance 
value of 3.71. The result of {f /g)uj = —0.06 is also found to be very close to the vector- 
meson dominance value of -0.12, but is not expected to hold exactly since the meson 
is not considered to be a pure ss system (we do not assume ideal mixing for the vector 
mesons). It is also interesting to note that the e coupling constant is much smaller than in 
the one-boson-exchange Nijm93 potential |T^, which is to be expected since in the present 
ESC model two-pion exchanges are explicitly included. The 14 free parameters are fitted to 
the 1993 Nijmegen representation of the hypersurface of the NN scattering data below 



^lab = 350 MeV ||28|, updated with the inclusion of new data which have been published 
since then. A comparison between this updated partial-wave analysis and the Nijm93 and 
ESC potentials will be given in the companion paper [0]. 

In Figs. ^ and ^ we compare the contributions from one-boson exchange, pion-pion 
exchange, and the other pion-meson exchanges; all for isospin 1 = and / = 1, respectively. 
The OBE contributions are the standard OBE potentials, substituting the coupling constants 
and cutoff masses as given in Table |IV|. Clearly, these OBE contributions are different from 
what one would obtain in a potential model containing only OBE contributions: the presence 
of the two-meson contributions substantially modifies the short-range behavior which has 
to be compensated for. Hence, the OBE contributions of the present ESC model on their 
own will certainly not fit the data. All potentials are seen to level off towards moderate 
values at the origin. This is due to the Gaussian form factors. (Note that the Gaussian 
form factors ensure that all potentials are finite at the origin to all orders in k^/M^ and 
that the tensor components always vanish at the origin; this in contrast to a monopole or 
dipole form factor.) The pion-meson contributions dominate the inner region; especially the 
spin-spin potential, which is mainly due to pion- vector exchange (see Figs. ^ and |^). For 
7 = 1 the one-boson and pion-pion central contributions cancel each other to a large extent, 
whereas for / = they enhance each other. The OBE contributions clearly dominate the 
outer region (r ^ 1 fm). 

In Figs. 1^ and ^, we compare the contributions to the pion-meson potentials from the 
different types of mesons: pseudoscalar (tt, rj, rj'), vector (p, u, 0), scalar (ao, e, /o), and 
Pomeron. The pion-scalar and pion-Pomeron potentials do not have a central component, 
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while the central component of the ttt] and nr]' potentials is very small, as can be inferred 
from the comparison of the pseudoscalar curve in Figs. |] and ^ and the pion-pion curve in 
Figs. ^ and 0, respectively. A similar comparison for the spin-spin and tensor components 
shows that, at short distances, the irr] and in]' potentials oppose the vrvr potentials for 
/ = 0, resulting in a largely reduced pion-pseudoscalar contribution. For 1 = 1 they slightly 
enhance the hit potentials. Looking into even more detail (not shown in the figures), we 
find that the ttt] and irrj' potentials are the least important. This is due to the small NNr) 
and NNr]' coupling constants and the high mass of the r]'. 

The large spin-spin pion-vector potential can be traced to the rather large uj coupling 
constant and cutoff mass. At the origin, the tcuj component reaches a maximum of 3122 MeV 
for / = 0, which is only partially reduced by the vrp component. For 1=1 the ttuj component 
reaches a minimum of -1040 MeV, which enhances the equally large up component. The ttcJ) 
component contributes much less, as is to be expected due to the large 0-meson mass. At 
present, it is not yet clear to us why the model apparently requires such a large spin-spin 
component. 

Comparing our results to those of Holinde and Machleidt (HM) |T^ we also find a large 
cancelation between the tttt and irp potentials. The cancelation between the ne and vru; 
potentials, however, is much less pronounced. One reason for this is the enormous difference 
in coupling constants: we have gl/^'n = 0.34 and g'^/4:TT = 11, whereas HM use 6 and 23, 
respectively. As a matter of fact, in our model the contribution of the ire potential is almost 
negligible. Another reason for the only partial cancelation is that we include the scalar /o 
and vector mesons, whereas HM do not. 

The 1/M contributions to the potential due to the nonadiabatic, pseudovector- vertex, 
and off-shell TMO 1/M corrections are shown in Fig. ^ Except for the off-shell TMO 
spin-spin part, all contributions are seen to vanish at the origin. Furthermore, there are 
large cancelations between the nonadiabatic and pseudovector-vertex corrections, and so 
the combined 1/M contribution to the pion- meson potential is indeed much smaller than 
the leading-order contribution, validating our procedure of expanding the potential as a 
series in 1/M. 

In Fig. 1^ we show the spin-orbit contributions of the pion-meson potentials in comparison 
with the spin-orbit OBE potential. The spin-orbit potentials of the pion-scalar and pion- 
vector exchanges are due to the 1/M^ terms in Eqs. ( [4.13| ) and ( [4.26| ). The spin-orbit 
pion-pseudoscalar potential, on the other hand, is mainly due to the 1/M pseudovector- 
vertex contributions, and hence is much larger. For / = 1, the pion-pseudoscalar spin-orbit 
potential largely cancels the OBE spin-orbit potential in the inner region, whereas for / = 
they enhance each other. The resulting spin-orbit part is reduced by the pion-vector spin- 
orbit contribution. (Incidentally, the large uj coupling constant and cutoff mass might be 
a necessary condition to arrive at a sufficiently moderate 1 = spin-orbit potential.) On 
the other hand, the 1/M nonadiabatic and off-shell TMO spin-orbit contributions, though 
much smaller than the pseudovector-vertex contribution, peak near 0.4 fm and, therefore, 
cause a substantially slower fall-off of the spin-orbit component. Because the spin-orbit 
operator L-S is proportional to the angular momentum L, this means that this contribution 
becomes more and more important for the higher partial waves, i.e., at higher energies. We 
indeed find that the quality of the fit rapidly becomes worse beyond ~300 MeV. It has 
been argued PDIISO] that these spin-orbit contributions will be (partially) cancelled by the 
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inclusion of higher-order (three-meson- and four-meson-exchange) diagrams, and so perhaps 
they should not be included at the present level of our expansion. Alternatively, one can 
argue that the nonadiabatic expansion in principle breaks down at the pion-production 
threshold (~280 MeV), and so the energy range should not be extended to 350 MeV anyway. 
However, the inclusion of the nonadiabatic central, spin-spin, and tensor contributions does 
provide a considerable improvement beyond ~280 MeV, apparently contradicting this line 
of argument. Currently, we are still investigating this matter. 

Finally, we note that if we choose to take the pseudoscalar coupling ( |2.13|) instead of 
the pseudovector coupling ( ^.12a| ), there are some small differences in the pion-scalar and 
pion- vector potentials; i.e., one of the terms changes sign, see Table |1[ A much larger 

effect is due to the fact that in case of the pseudoscalar coupling the u/M pseudovector- 
vertex corrections of Sec. |V B| are all absent. Also, there is a sign change in the off-shell 
TMO corrections, which means that in the pion-pseudoscalar potential the nonadiabatic and 
off-shell TMO spin-orbit contributions exactly cancel (they have the same r dependence). 
Hence, the apparent problem with the spin-orbit contribution as discussed above does not 
occur. 
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APPENDIX A: DIFFERENTIATION DICTIONARY 



In this Appendix we give a dictionary for the evaluation of the differentiations in 
Eqs. ([4.71), ( [4.131 ), ( [4.26[ ), and the expressions in Sec. |V]. The procedure is described in 
more detail in Appendix B of [Q. In the following it is to be understood that the functions 
on the right-hand side are functions of r and that the prime denotes differentiation with 
respect to r. 

1. For the pseudoscalar operators Op^: 

lim (Vi-V2)'F(ri)G(r2) = ^F'G' + F"G\ 



lim ((Ti-VixV2)(<T2-VixV2)F(ri)G(r2) = - 

ri^r2 3 

1 

~3 



Lf'G' + -F'G" + -F"G' 
F" - -F'\ -G' + -F' ( G" - -G' 



(Ti-(T2) 

Sl2- 



2. For the scalar operators O^^ and the vector operators O^*^, O^j^, O'^^lm^ ^^^d 



' m,m' 



lim (o-rVi)(o-2-Vi)F(ri)G(r2) = - 

ri^r2 3 

lim (<TrVi)((T2-Vi)V^F(ri)G'(r2) = 



F" + ^F') G{a,-(T2) + (f" - ^F' ) G Si2 
F" + V) (g" + Ig') (.,.2) 
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+ [F"- -F' ]{G" + -G' ] S_ 



r 



r 



'12 



hm^(o-i-Vi)((T2-Vi)(Vi-V2)F(ri)G'(r2) = - 



(o-i ■ Vi) (0-2 • V2) V2F(ri)G'(r2) 
^Im (o-i ■ V2) (0-2 ■ Vi) V?F(ri)G'(r2) 
^Im ((Ti ■ Vi) (0-2 ■ Vi) V?F(ri)G'(r2) 



+ F 



2 2 
2 2 
4 



G"[(o-ro-2) + 5i2] 



F"" + -F'")GicT^■(T2 



F"" + If'" - ^F" + ) G S 



12 



lim [q-Vi X V2] {(Ti + (T2)-ViF{n)G{r2) = -F'G' L-S. 



ri^r2 



3. For the nonadiabatic and pseudovector- vertex corrections: 



6 



ri^r2 



1 



lim (Vr V2)=^F(ri)G'(r2) = -iF"--F')i G" - -G' ) + F"'G 



1 



^lim^( Vi ■ V2) (cTi ■ Vi X V2) (0-2 ■ Vi X V2)F(ri)G(r2) 



2 
1 

+3 



i;(-F'-F" + rF'" ) ( -G" - G" + rG 



ty'2i \ 



F"'G" 
1 



Icri-cr2 



\(-F'- 2F" + '-F'" U -G" - 2G" + ^G'" ) - ^F"'G 



S 



12, 



^lim^(Vi- V2)(o-i + tT2)-(Vi X V2) [q-(Vi - V2)]F(ri)G(r2) 



('g'"--g"')l-s, 

r J \ r J 



Jim (^Tr Vi)(o-2- V2)F(ri)G(r2) = ^FG" [(tTro-2) + S12] , 



ri^r2 



^lim^(<Ti- V2)(o-2- Vi)F(ri)G(r2) = ^F'G' [{(T^-a2) + ^12] , 

^lim^(cri + cr2)-q x (Vi + V2) (Vi ■ V2)F(ri)G(r2) = [F"G' + F'G"] L-S, 

^lim^(cri + cr2)-(Vi x V2) [q-(Vi - V2)]F(ri)G(r2) = -^F'G'L-S. 



APPENDIX B: 
COORDINATE-SPACE POTENTIALS 



In this appendix we give the exphcit expressions for the coordinate-space two-meson- 
exchange potentials. We also include the result for the two-pion-exchange potential from 
Ref. for reasons of completeness. In order to keep the expressions as compact as possible, 
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we restrict ourselves to the leading-order terms (no contributions), and define and 

operations as 

(F0G).(r) = i(F" + V) G(r), 
(F0G)T(r) = i(F"- V) G(r), 



(F G)c;(r) 
(F0G).(r) 
(F G')T(r) 



+ —F'G' 



^F'G' + + -F'G" 

^F'G' - -F"G' - -F'G" 



(r). 



We then find the following [FQ,(A,r) is given by Eq. ( [4.5|) ]. 
1. Pseudoscalar exchanges: 



TTTT 



n Jo 



X 2(TrT2) 



m - {F. F.; 



c 



+3 



+3 



m - {F^ F^)^ 



CTi-CTal 



5'l2 



\/(7rr7) 



V / \ / vr JO 



(^2 ® n). - {F, ® F,x 



l(Ti-(T2 



+ 



I-Dt ~ (F^ 



vJt 



S 



12 



V{TTr]') = ¥{717]) with rj rj'. 



(Bl) 



(B2) 
(B3) 



2. Scalar exchanges: 



^(vrao) = 4^, 



2 / /aWtt 



,2 /•°°(iA 

T1-T2 - / 



[Cri-(T2 



V{7ie) 
3. Pomeron exchange: 



ViTTfo) 



2 ;t 
0. 



{F^ 



5 



12 f , 



(B4) 
(B5) 
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V{TrP) = 0. 



(B6) 



4. Vector exchanges: 



V{np) = -Agl 



'I 



NNn 



,2 f^dX 

TrT2 - / 



X 



+ 



S 



12 



l^(7rcu) = 1/(7r0) = 0. 



(B7) 
(B8) 
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TABLES 



TABLE I. Adiabatic approximation of the energy denominators Di{uji,uj2) for the planar and 
crossed diagrams. 



planar (//) 
crossed {X) 



1 



1 



1 


+ 


1 






1 


+ 


1 




0J2 



1 
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TABLE II. The momentum operators 0^*^(ki,k2) for the planar (//) and crossed {X) 
TT-meson-exchange diagrams. The subscripts PS and 5* refer to the vr-pseudoscalar and 7r-scalar 
operators, respectively. The subscripts e and m refer to the electric and magnetic parts of the 
TT- vector operators, and k = fv/dv- The ^ sign refers to pseudovector or pseudoscalar coupling 

for the pioii, respect iv(;ly. 



Type 



OW(ki,k2) 



Oj//) 2(ki •k2)2 - 20-1 • (ki X k2)o-2 • (ki X k2) 
2(ki •k2)2 + 2(Ti ■ (ki X k2)(T2 • (ki X k2) 



o 



ill) 



o 



2 + 



ki+krk2 q2 + ik2 



2M2 



4M2 



M2 



criki(T2ki 



{kj ^ ki-k2)(<Ti-ki<T2-k2 + o-i-k20-2-ki) + 2iq-(ki xk2)((Ti + cr2)-ki 



2 + 



ki+ki-k2 q2 + ik2 



2M2 



4M2 



M2 



a-ikio-2ki 



(ki ^ ki-k2)((Ti-ki<T2-k2 + <Ti-k2(T2-ki) + 2zq-(ki xk2)(<Ti + (T2)-ki 



-2-3 



q2 + ik2 ^ k2 + 4krk2 + k2 



M2 

.2 , ll.2 



+ 



2M2 



2M2 



^ ^ q^ + jk^ ^ k2 + 4ki-k2 + k2 ^ k? 



M2 



2M2 



2M2 



o-ikio-2ki =F 

<Ti ki(T2 ki T 



kik2 

4M2 

kik2 

4M2 



(cri-kicr2-k2 + (7i-k2cr2-ki) 

(<Ti-ki(T2-k2 + (Ti-k2(T2-ki) 



(//) 



(1 + k) [-27;q-(ki xk2)(cri + <T2)-ki + 2ki-k2cri-kicr2-ki - kf (cri-kicr2-k2 + <ti •k2cr2-ki)] 
(1 + k) [-2iq-(ki xk2)(o-i + cr2)-ki - 2ki-k20-i-kicr2-ki + k2(cri-kicr2-k2 + cri •k20-2-ki)] 



OiilL (1 + 2k) [-iq-(ki Xk2)(<Ti + tT2)-ki + (k| + ki-k2)<Tl-ki(T2-ki 



k2(cri-kicr2-k2 + cri •k2cr2 -ki) 



Oifn (1 + 2«) [-iq- (ki X k2)(ai + 0-2) -ki + (k2 + ki •k2) o-ikio-2ki 



-ikf (o-i-ki(T2-k2 + cri-k2(T2-ki) 



Omfii (1 + k)2 [kfk^ - (ki-k2)^ + k|o-i-kio-2-ki - kfk|o-i-(T2 + o-i-(ki xk2)c7-2-(ki xk2 



O^i (1 + [kf k| - (ki •k2)2 - k^o-i •kio-2 -ki + kf k^tri •0-2 - <ti • (ki x k2)o-2 • (ki x k2)] 
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TABLE III. Energy denominators dI^\uji,uj2) ^ov the nonadiabatic corrections to the planar 
and crossed diagrams. 



planar (//) 
crossed {X) 



+ 



1 



1 



+ 



1 1 



TABLE IV. 


Meson parameters 


employed in 


the potentials shown in 


Figs. 1 to 1 


. Coupling 


constants are at = 0. 










meson 


mass (MeV) 


g/y/4TT 




f/g 


A (MeV) 




138.04 




0.269 




853.9 




547.45 




0.124 




950.0 


rj' 


957.75 




0.124 




950.0 


P 


768.10 


0.622 


2.306 


3.71 


886.2 




781.95 


3.336 


-0.206 


-0.06 


1013.1 




1019.41 


-1.202 


0.078 


-0.07 


1013.1 


ao 


982.70 


1.847 






679.4 


e 


760.00 


0.582 






759.0 


fo 


974.10 


-2.372 






759.0 


Pomeron 


309.10 


1.884 
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FIGURES 



FIG. 1. BW two-meson-exchange graphs: (a) planar and (b)-(d) crossed box. The dotted hne 
with momentum ki refers to the pion and the dashed hne with momentum k2 refers to one of the 
other (vector, scalar, or pseudoscalar) mesons. To these we have to add the "mirror" graphs, and 
the graphs where we interchange the two meson lines. 

FIG. 2. Planar-box TMO two-meson-exchange graphs. Same notation as in Fig. I. To these 
we have to add the "mirror" graphs, and the graphs where we interchange the two meson lines. 

FIG. 3. Central, spin-spin, and tensor potentials for / = 0. Shown are the one-boson exchange 
(OBE), the pion-pion exchange, and the other pion- meson exchanges. All potentials are for (a) 
r < 1 fm and (b) 1 < r < 2 fm. The spin-spin pion-meson potential has its maximum at 4874 
MeV. 

FIG. 4. Same as Fig. |3|, but for 1 = 1. The spin-spin pion-meson potential has its minimum at 
-2969 MeV. 

FIG. 5. Central, spin-spin, and tensor potentials for 7 = 0. Shown are the pion-pseudoscalar, 
pion-vector, pion-scalar, and pion-Pomeron exchanges. All potentials are for (a) r < 1 fm and (b) 
1 < r < 2 fm. 

FIG. 6. Same as Fig. |, but for / = 1. 

FIG. 7. Nonadiabatic, pseudovector-vertex, and purely off-shell 1/M contributions to the 
pion-meson potentials for (a) 1 = and (b) 1 = 1. 

FIG. 8. Spin-orbit contributions of the pion-meson potentials for both 7 = and / = 1 in 
comparison with the one-boson-exchange contribution. 
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Figure 1 , paper I ( Rijken, Stoks ) 
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Figure 



2, paper I ( Rijken, Stoks ) 
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Figure 3, paper I (Rijken, Stoks) 




Figure 3, paper I (continued) 




Figure 4, paper I (continued) 




Figure 5, paper I (Rijken, Stoks) 




Figure 5, paper I (continued) 




Figure 6, paper I (Rijken, Stoks) 
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Figure 6, paper I (continued) 




Figure 7, paper I (Rijken, Stoks) 




Figure 7, paper I (continued) 
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Figure 8, paper I (Rijken, Stoks) 



